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Abstract. This paper describes the comparison of the vari- 
ability of total column ozone inferred from the three inde- 
pendent multi-year data records, namely, (i) Solar Backscat- 
ter Ultraviolet Instrument (SBUV) v8.6 profile total ozone, 
(ii) GTO (GOME-type total ozone), and (iii) ground-based 
total ozone data records covering the 16-year overlap period 
(March 1996 through June 2011). Analyses are conducted 
based on area-weighted zonal means for 0-30° S, 0-30° N, 
50-30° S, and 30-60° N. 

It has been found that, on average, the differences in 
monthly zonal mean total ozone vary between —0.3 and 
0.8 % and are well within 1 %. 

For GTO minus SBUV, the standard deviations and ranges 
(maximum minus minimum) of the differences regarding 
monthly zonal mean total ozone vary between 0.6-0. 7 % and 
2. 8-3. 8 % respectively, depending on the latitude band. The 
corresponding standard deviations and ranges regarding the 
differences in monthly zonal mean anomalies show values 
between 0.4-0. 6% and 2. 2-3. 5%. The standard deviations 
and ranges of the differences ground-based minus SBUV re- 
garding both monthly zonal means and anomalies are larger 
by a factor of 1.4-2. 9 in comparison to GTO minus SBUV. 

The ground-based zonal means demonstrate larger scatter- 
ing of monthly data compared to satellite-based records. The 
differences in the scattering are significantly reduced if sea- 
sonal zonal averages are analyzed. 


The trends of the differences GTO minus SBUV and 
ground-based minus SBUV are found to vary between —0.04 
and 0.1 %yr -1 (—0.1 and 0.3DUyr _1 ). These negligibly 
small trends have provided strong evidence that there are 
no significant time-dependent differences among these multi- 
year total ozone data records. 

Analyses of the annual deviations from pre-1980 level in- 
dicate that, for the 15-year period of 1996 to 2010, all three 
data records show a gradual increase at 30-60° N from —5 % 
in 1996 to -2 % in 2010. In contrast, at 50-30° S and 30° S- 
30° N there has been a levelling off in the 15 years after 1996. 
The deviations inferred from GTO and SBUV show agree- 
ment within 1 %, but a slight increase has been found in the 
differences during the period 1996-2010. 


1 Introduction 

Concern for changes in the ozone layer due to human ac- 

tivity is an important subject for the scientific community, 

the general public and governments. Accurate long-term data 

records of total column ozone and vertical profiles of ozone 

are required for the scientific assessment of ozone depletion 

(WMO, 201 1). 


Published by Copernicus Publications on behalf of the European Geosciences Union. 


1682 


E. W. Chiou et al.: Comparison of profile total ozone from SBUV (v8.6) 


In response to the observed ozone loss, countries around 
the world adopted the Montreal Protocol and subsequent 
amendments calling for limitations on production and use of 
ozone-depleting substances (UNEP, 2006; UNFCCC, 1998). 
Besides the impact of ozone depleting substances, natural 
fluctuations such as the 11 -year solar cycle, the equatorial 
quasi-biennial oscillation of the lower stratospheric zonal 
wind (QBO), and volcanic eruptions also significantly affect 
the thickness of the ozone layer. Moreover, climate change 
due to increase in greenhouse gas concentration will influ- 
ence stratospheric dynamics and chemistry and therefore the 
ozone layer. Many investigations have been conducted for 
monitoring and detection of global ozone trends and behav- 
ior using a variety of ground-based and satellite instruments 
and their comparisons. Recent studies based on long-term 
ozone data records and model simulations have significantly 
improved our understanding in the roles of various dynam- 
ical and chemical processes governing the ozone variations 
(e.g., Yang, 2006; Stolarski and Frith, 2006; WMO, 2011). 
However, many detail characteristics of the expected ozone 
recovery such as the beginning of the recovery and the timing 
of the recovery are still unclear. One of the major difficulties 
in assessing long-term global ozone variations is data inho- 
mogeneity. Changes in operational satellites, revision of re- 
trieval algorithms, recalibration of ground-based instruments 
or interruptions in observation periods result in data sets that 
have systematic errors that change with time. 

The purpose of this study is to conduct an investiga- 
tion of the consistency in the variability of global and 
zonal total ozone inferred from three independent multi- 
year data records, namely, (1) the recently released So- 
lar Backscatter Ultraviolet Instrument (SBUV) v8.6 pro- 
file total ozone data record, (2) the new European GOME- 
type total ozone record (GTO), and (3) ground-based to- 
tal ozone data record based on Dobson and Brewer spec- 
trometer and filter ozonometer observations available from 
the World Ozone and Ultraviolet Radiation Data Cen- 
tre (WOUDC) (http://www.msc-smc-smc.ec.gc.ca/woudc/). 
The two satellite -based merged ozone data records (SBUV 
v8.6 and GTO) are the most recent versions released during 
2012-2013. They represent the results of two independent 
attempts by US and European scientists to adopt an opti- 
mal technique to merge ozone measurements from a series 
of satellite instruments aiming to construct a homogeneous 
self-consistent and calibrated long-term data record. 

The analyses are based on the 16-year overlap period of 
March 1996 through lune 2011. Our major goal is to quan- 
titatively evaluate the long-term stability of the global and 
zonal mean total column ozone from these three multi-year 
data records. Discussions will be focused on the consistency 
regarding monthly zonal mean total ozone and the monthly 
zonal mean anomalies. The results obtained herein should 
enable us to enhance our understanding of the accuracy of 
these multi-year data records. 


Detailed information of the three data records is pro- 
vided in Sect. 2. Results of comparisons based on 5-degree 
monthly zonal means between 60° S and 60° N are presented 
in Sect. 3.1. Further analysis using area-weighted monthly 
zonal means for 0-30° S, 0-30° N, 50-30° S, and 30-60° N 
is discussed in Sects. 3.2 and 3.3. Changes that will occur 
when switching to analysis using seasonal zonal means are 
reported in Sect. 3.4. Trends of the differences are discussed 
in Sect. 3.5. Investigation of total ozone deviations from pre- 
1980 levels is reported for 30-60° N, 50-30° S, and 30° S- 
30° N in Sect. 3.6. Concluding remarks from our study are 
summarized in Sect. 4. 


2 Description of the three ozone data records 

2.1 SBUV v8.6 profile total ozone data record 

NASA and NOAA have been measuring ozone from space 
since 1970. The previously existing merged ozone data set 
provided by NASA combines the TOMS data (Nimbus 7 and 
Earth Probe) and SBUV-SBUV/2 data (Nimbus 7, NOAA 9, 
1 1, 14, 16). Studies were made by Stolarski and Frith (2006) 
using the merged ozone data set to search for evidence of 
ozone recovery in response to the observed levelling off of 
chlorine compounds in the stratosphere. Before v8.6, the phi- 
losophy for producing the merged data set was to take the in- 
dividual data sets and combine them by making simple off- 
set corrections in ozone based on overlap periods or com- 
parisons with other data sets. The correction offset is deter- 
mined as the average difference in the 50° S to 50° N zone. 
No time dependence is applied to an individual data set. All 
the monthly zonal means for an instrument are adjusted by a 
single offset. There is no latitudinal dependence applied. All 
data sets available for each month are then averaged together 
to produce the final merged ozone time series. 

However, for v8.6 no global offsets or any offsets were 
applied to profile total ozone records from individual SBUV 
instruments before combining them. The inter-calibration is 
accomplished within the v8.6 algorithm (Frith et al., 2014). 
The SBUV v8.6 profile total ozone data record incorporated 
the measurements from eight backscatter ultraviolet instru- 
ments (BUV on Nimbus 4, SBUV on Nimbus 7, and a se- 
ries of SBUV/2 instruments on NOAA satellites). The cov- 
erage periods of each instrument used to create the merged 
ozone data set are shown in Table 1. As described in Frith 
et al. (2014), pairs of the SBUV measurements records are 
combined using a simple average during periods when data 
from more than one instrument are present. The only ex- 
ception is the exclusion outside the ECT (Equator crossing 
time) range of 8 a.m. to 4 p.m. and nearly all NOAA-9 data. 
Discussion of SBUV v8.6 algorithm is presented by Bhartia 
et al. (2013). An overview of the version 8.6 SBUV ozone 
data record is discussed by McPeters et al. (2013). Major im- 
provement has been achieved by radiance adjustment made 
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Table 1. Period of coverage by each instrument used in constructing the SBUV merged profile total ozone data record. 


Satellite instrument Period covered Remarks 


Nimbus 4 BUV 
Nimbus 7 SBUV 
NO A A 1 1 SBUV/2 
NO A A 1 1 SBUV/2 
NO A A 14 SBUV/2 
NO A A 14 SBUV/2 
NO A A 16 SBUV/2 
NO A A 16 SBUV/2 
NO A A 17 SBUV/2 
NO A A 18 SBUV/2 


May 1970- Apr 1976 
Nov 1978-May 1990 
Jan 1989-Mar 1995 
Oct 1997-Mar 2001 
Mar 1995- Apr 2000 
Jul 2004-Sep 2006 
Oct 2000-Jun 2007 
Dec 2011-Dec 2011 
Aug 2002-Sep 2011 
Jul 2005-Dec 2011 


(All data) 

(All data) 

(Enter terminator orbit; 18:00ECT) 
(08:00 ECT; end of record) 

(16:00 ECT; nearing terminator) 
(08:00 ECT; end of record) 

(16:00 ECT; nearing terminator) 
(08:00 ECT) 

(08:00 ECT; nearing terminator) 


* ECT: Equator crossing time. 


for each instrument to maintain a consistent calibration (De- 
Land et al., 2012). This new merged ozone data record covers 
the period from 1970 to 2011. Two other important changes 
in the processing of v8.6 are that the ozone cross sections of 
Brion, Daumont and Malicet have been used and that a cloud 
climatology derived from the Aura/OMI (Ozone Monitoring 
Instrument) cloud-height retrievals has been used (Bhartia et 
al., 2013). 

Results of the study for the effect of applying the new cross 
sections to SBUV retrievals indicate that the ozone profiles 
are lower in the upper stratosphere and higher in the lower 
stratosphere and troposphere. As a result, the average differ- 
ences in total ozone are small, and close to zero in the tropics 
(McPeters and Labow, 2010). 

Recent studies of the differences in total ozone between 
v8.6 and v8.0 (Merged Ozone Data Sets) MODs indicate 
that v8.6 is 0.5 to 1 % lower than v8.0 at all times except 
in the mid-1990s, when v8.6 is 1 % higher. The difference in 
mid- 1990 is related to the use of NOAA-14 in the new ver- 
sion rather than the NOAA-9 used in the v8.0 MOD product. 
The deviation in mid- 1982 results from different treatment of 
the data after the eruption of El Chichon in late March 1982 
(Frith et al., 2014). 

2.2 GTO ozone data record 

The European satellite-borne sensors GOME/ERS-2 (1995- 
2011), S Cl AM ACHY /EN VIS AT (2002-2012), and GOME- 
2/METOP-A (2007-present) provide global total ozone mea- 
surements for the last 17 yr. A summary of the instrument 
properties and viewing geometries is given in Table 2. 

The GTO merged ozone data record combines those mea- 
surements, and a continuous and homogeneous monthly 
mean time series is generated (Loyola et al., 2009; Loy- 
ola and Coldewey-Egbers, 2012). The first GTO version 
was created using products obtained with the GOME Data 
Processor (GDP) version 4.x algorithm (Van Roozendael et 
al., 2006; Lerot et al., 2009; Loyola et al., 2011), which 
is based on the differential optical absorption spectroscopy 


approach. Geophysical validation shows that GDP 4.x total 
ozone has an accuracy at the percentage level compared with 
ground-based instruments (Loyola et al., 2011; Koukouli et 
al., 2012). The resulting GTO data record was used for the 
WMO ozone assessment report 2010 (WMO, 2011) and for 
ozone studies (Dameris and Loyola, 2012). 

In this study we use the most recent version of the GTO 
data record that has been developed within the framework 
of the European Space Agency’s Climate Change Initiative 
(ESA-CCI). It incorporates the ozone data products retrieved 
using the newly developed GOME Direct Fitting algorithm 
GODFIT (Lerot et al., 2010, 2014; Van Roozendael et al., 
2012) and covers the period from March 1996 to June 2011. 

The GTO merging approach accounts for the generally 
small remaining differences among the individual instru- 
ments, which mainly depend on latitude, season, and time. 
Due to excellent long-term stability, the GOME measure- 
ments are used as a transfer standard, whereas SCIAMACHY 
and GOME-2 data are adjusted accordingly in periods of in- 
strument overlap (Loyola et al., 2009). The adjustments com- 
prise two parts: a basic latitudinal correction for each month 
of the year averaged over all years and a time-dependent off- 
set for each individual month. Within the ESA-CCI the merg- 
ing approach has been further improved; the adjustment is 
now based on daily averages (instead of monthly) in order to 
reduce differences among the individual instruments due to 
differences in spatial and temporal sampling (see Table 2). 
The daily averages now incorporate all available measure- 
ments from one day, which increases the number of observa- 
tions at high latitudes where the orbits overlap. Monthly av- 
erages are computed taking into account latitudinal and sea- 
sonal constraints in order to provide representative mean val- 
ues. Finally, GOME, adjusted SCIAMACHY and adjusted 
GOME-2 data are merged. 

2.3 Ground-based ozone data record 

The World Ozone and Ultraviolet Radiation Data Centre 
(WOUDC) in Toronto began collecting and publishing ozone 
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Table 2. Characteristics of the three European satellite instruments. 


Parameter 

GOME 

SCIAMACHY 

GOME-2 

Data availability 

Jul 1995-Jun 2011* 

Aug 2002-Mar 2012 

Jan 2007-today 

Spectral coverage 

240-790 nm 

240-2380 nm 

240-790 nm 

Spectral resolution 

0. 2-0.4 nm 

0.2-1. 5 nm 

0. 2-0.4 nm 

Ground pixel size 

320 x 40 km 2 

60 x 30 km 2 

40 x 80 km 2 

Swath width 

960 km 

960 km 

1920 km 

Equator crossing 

10:30 LT 

10:00 LT 

09:30 LT 

Global coverage 

3 days 

6 days 

almost daily 


* No global coverage since June 2003. 


Geographical Distributions of WOUDC stations 



Dobson Stations 
Brewer Stations 
Filter Instrument Stations 

Figure 1. Geographical distributions of WOUDC ground stations. 

data in 1961 and remains the main source of ground-based 
ozone data for researchers. Global distributions of ground- 
based stations used in this study are depicted in Fig. 1. 
Large longitudinal inhomogeneity and limited spatial cov- 
erage make it impossible to estimate zonal and global to- 
tal ozone values from station values directly. However, if 
an ozone “climatology” (i.e., long-term mean for each point 
of the globe for each day of the year) estimated from 
satellite data is used with ground-based measurements of 
ozone deviations from that climatology at the stations, then 
long-term zonal and global ozone variations can be estimated 


using ground-based data (Bojkov and Fioletov, 1995). The 
present climatology used is based on TOMS N7 version 8 
for the 1978-1987 period everywhere except for the Antarc- 
tic, where 1978-1982 data only were used. The climatology 
was last updated in 2004. In summary, the method measures 
ozone deviations from the climatology at the stations, then 
calculates the zonal deviations, and finally the zonal mean 
ozone is determined by adding the zonal mean climatology 
to the zonal means of the deviations. The results give a con- 
tinuous uninterrupted global total ozone data record that is 
fairly independent of other data sources. 

However, the absence of data over vast regions (e.g., 
oceans) and sensitivity to individual instrument errors is an 
important factor, particularly in the tropical region and the 
Southern Hemisphere where the number of stations is very 
limited. Comparisons of global and zonal ozone variations 
from ground-based and satellite measurements for the pe- 
riod 1964-2000 were presented by Fioletov et al. (2002). 
The ground-based multi-year total ozone record used in this 
study is the recently updated one extending through Novem- 
ber 2012. 

3 Results and discussions 

Even though the most recently released SBUV (v8.6) ozone 
profile and ground-based total ozone data records cover a 
43-year period (1970-2012) and a 49-year period (1964- 
2012), respectively, our studies focus on the 16-year period 
(March 1996 to June 2011) for which these two data records 
overlap with GTO. For high-latitude regions, the number of 
ground stations is very limited to represent the characteris- 
tics of the latitudinal zones. It is also known that the available 
number of SBUV measurements is drastically reduced in the 
high-latitude regions in both hemispheres. For example, as 
shown in Fig. 2, among the 184 months (from March 1996 
to June 2011) selected for this study, more than 20% of 
the months are missing data poleward of 60° S and 65° N in 
SBUV (v8.6) data records. Thus, our analyses and compar- 
isons will only be based on the 24 5-degree monthly zonal 
mean time series (covering 60° S to 60° N) from these three 
multi-year data records. 
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Missing Data in High Latitudes 



Figure 2. Missing data in high latitudes for SBUV v8.6. 


3.1 The 5-degree monthly zonal means between 60° S 
and 60° N 

Figure 3a-c show the monthly mean zonal mean total column 
ozone as a function of latitude and time given by SBUV v8.6, 
GTO and ground-based data records, respectively. Black re- 
gions in Fig. 3a represent missing data in SBUV. The fea- 
tures revealed in these figures clearly indicate that the three 
data records exhibit almost identical patterns of temporal and 
latitudinal variations for the entire 16-year period. 

The differences (in %) between each pair of data records 
are illustrated in Fig. 4a-c. Black regions in Fig. 4a and b 
represent missing data in SBUV. The agreement between the 
two satellite-based ozone data records is significantly better 
than the agreement between each satellite data record and 
the ground-based record. Statistical parameters for the per- 
centage differences revealed in these figures are listed in Ta- 
ble 3a-c. 

After a more detailed investigation of the three data 
records, it has been found that the monthly zonal mean at 
60-55° S from ground-based record has an obvious outlier 
appearing at one particular month - November 2009. Thus, 
we have intentionally excluded that particular data point in 
deriving the statistical parameters listed in Table 3b and c. 
The entries with an asterisk listed in the first row in Table 3b 
and c would increase to 27.6 and 28.3 % respectively without 
excluding this data point. 

Major findings from Table 3a-c can be summarized as fol- 
lows: 


(a) SBUV Data Record 



(b) GTO (GOME Type total ozone version 2) 



1997 1999 2001 2003 2005 2007 2009 2011 

16-year period (MAR 1996-JUN 2011) 


total ozone (DU) 



229 283 337 391 445 

(c) Ground-based Data Record 


n ' ' i rr 



JAN JAN JAN JAN JAN JAN JAN JAN 

1997 1999 2001 2003 2005 2007 2009 2011 


16-year period (MAR 1996-JUN 201 1) 


total ozone (DU) 



229 283 337 391 445 


Figure 3. Monthly mean zonal mean total column ozone (5- 
degree zones): (a) SBUV, (b) GTO (GOME-type total ozone), and 
(c) ground-based data record. (Black regions represent missing data 
in SBUV v8.6.) 


1. For both GTO versus SBUV and satellite-based record 
versus ground-based record, the standard deviations and 
the ranges (maximum minus minimum, representing 
peak-to-peak variations) of the differences are simi- 
lar for 0-30° S, 0-30° N, and 30-60° N and are much 
smaller than those at 30-60° S. 

2. The standard deviations of the differences between the 
two satellite -based records vary between 0.7 and 1.0 %. 


The corresponding standard deviations for satellite- 
based record versus ground-based record are signifi- 
cantly larger, ranging between 1.4 and 2.1 %. 

3. The ranges of the differences GTO minus SBUV vary 
between 3.9 and 9.9%. The corresponding ranges of 
the differences for ground-based record versus satellite- 
based record are larger by a factor of 2 to 3. 
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(a) GTO minus SBUV 



16-year period (MAR 1996-JUN 2011) 

Differences (%) 



(b) Ground-based record minus SBUV 



JAN JAN JAN JAN JAN JAN JAN JAN 

1997 1999 2001 2003 2005 2007 2009 2011 


16-year period (MAR 1996-JUN 2011) 

Differences (%) 



-4 -2 -1 1 2 4 

(c) GTO minus Ground-based record 



JAN JAN JAN JAN JAN JAN JAN JAN 

1997 1999 2001 2003 2005 2007 2009 2011 


16-year period (MAR 1996-JUN 201 1) 

Differences (%) 



-4 -2 -1 1 2 4 


Figure 4. Differences in monthly zonal means (5-degree zones): 
(a) GTO minus SBUV, (b) ground-based minus SBUV, and (c) GTO 
minus ground-based. (Black regions represent missing data in 
SBUV v8.6.) 


Table 3. Summary of the differences in monthly zonal mean total 
ozone. 


(a) GTO minus SBUV 

Latitudinal zones 



Differences in percent 


Mean 

Standard deviations Range (Max-Min) 

(30-60° S) 

0.55 


0.98 

9.87 

(0-30° S) 

0.74 


0.75 

5.65 

(0-30° N) 

0.43 


0.70 

3.96 

(30-60° N) 

0.18 


0.78 

6.47 

(b) Ground-based minus SBUV 

Latitudinal zones - 



Differences in percent 


Mean 

Standard deviations Range (Max-Min) 

(30-60° S) 

-0.61 


2.07 

19.17* 

(0-30° S) 

0.56 


1.37 

11.32 

(0-30° N) 

0.10 


1.67 

12.71 

(30-60° N) 

0.67 


1.56 

12.04 

(c) GTO minus ground-based 

Latitudinal zones 



Differences in percent 


Mean 

Standard deviations Range (Max-Min) 

(30-60° S) 

1.21 


2.13 

19.89* 

(0-30° S) 

0.19 


1.39 

12.72 

(0-30° N) 

0.33 


1.62 

12.42 

(30-60° N) 

0.11 


1.66 

12.18 


* One outlier data point from ground-based record (November 2009; 60-55° S) 
was excluded. 


3.2 Area-weighted monthly mean zonal means for 
broader latitudinal zones 

In order to examine the characteristics of the differences 
between pairs of data records for several wider latitudinal 
zones, we have computed the area-weighted monthly mean 
zonal mean time series for (i) 0-30° S, (ii) 0-30° N, (iii) SO- 
SO 0 S, and (iv) 30-60° N. For the Southern Hemisphere, area- 
weighted means are limited to equatorward of 50° S due to 
a large number of months (namely, 26 months) with miss- 
ing data in the SBUV v8.6 data record. The monthly mean 
zonal mean time series for these four broad latitudinal zones 
are plotted in Fig. 5. The corresponding differences, with re- 
spect to SBUV, are illustrated in Fig. 6. The statistical pa- 
rameters corresponding to these differences are summarized 
in Table 4. 

For the ground-based minus SBUV differences, there are 
two periods, namely, beginning of 2007 and end of 2009 
where the differences show larger positive deviations in the 
0-30° N band. Figure 6c and d also indicate that the differ- 
ences GTO minus SBUV show slight changes from negative 
difference in 1996 to positive differences in the second half 
of the period. An explanation of these features is not read- 
ily available and further investigation in this regard is worth- 
while once revised data records become available in the near 
future. 
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Table 4. Summary of the differences in monthly zonal mean total ozone (four broad latitudinal zones). 


Latitudinal GTO minus SBUV (in %) Ground-based minus SBUV (in %) 


zones 

Mean 

Standard 

deviation 

Range 

(Max-Min) 

Mean 

Standard 

deviation 

Range 

(Max-Min) 

(0-30° S) 

0.76 

0.63 

3.63 

0.56 

0.91 

5.14 

(0-30° N) 

0.43 

0.58 

2.83 

0.10 

1.12 

6.92 

(50-30° S) 

0.67 

0.66 

3.82 

-0.32 

1.23 

5.75 

(30-60° N) 

0.20 

0.63 

3.35 

0.09 

0.98 

6.31 


Monthly Zonal Mean total column ozone 





16-year period (MAR 1996 - JUN 2011) 


: GTO (GOME Type Total Ozone) 


SBUV(v8.6) 


□ □ □ Ground-based Records 


Figure 5. Monthly mean zonal mean total ozone (area-weighted 
zonal means for broader latitudinal zones). 


The entries of mean differences in Table 4 have led to the 
conclusion that, on average, the differences of monthly mean 
zonal mean total ozone between the two satellite records and 
between satellite record and ground-based record vary be- 
tween —0.3 and 0.8 % and are well within 1 %. The long- 
term stability of the differences during the 16-year period 
is revealed by the standard deviations and ranges listed in 
Table 4. The standard deviations for GTO minus SBUV are 
found to lie between 0.6 and 0.7 %, while ground-based mi- 
nus SBUV exhibits less stability with standard variations 
varying between 0.9 and 1.2 %. The ranges of the differences 
for GTO minus SBUV vary between 2.8 and 3.8 %. It is also 





16-year period (MAR 1996 - JUN 2011) 


GTO minus SBUV 


Ground-based minus SBUV 


Figure 6. Differences in monthly mean zonal mean total ozone. 


noticed that the corresponding ranges of the differences for 
ground-based minus SBUV are larger by a factor of 1 .4 to 
2.4, varying between 5.1 and 6.9 %. 

3.3 Monthly mean zonal mean anomaly 

The consistency among the data records in terms of interan- 
nual variations in total ozone can be examined through in- 
vestigation of the monthly mean zonal mean anomaly. For 
each of the three data records, the monthly mean zonal mean 
anomaly is calculated by simply subtracting the 16-year av- 
erage for the same calendar month from the monthly mean 
zonal mean total ozone. The times series of monthly mean 
zonal mean anomaly for the four broad zonal bands are plot- 
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1995 2000 2005 2010 

16-year period (MAR 1996 - JUN 2011) 


xxx GTO (GOME Type Total Ozone) SBUV(v8.6) 

□ □ □ Ground-based Records 

Figure 7. Monthly mean zonal mean anomaly. 


ted in Fig. 7. The corresponding differences, with respect 
to SBUV, are illustrated in Fig. 8. The statistical parame- 
ters corresponding to these differences are listed in Table 5. 
The standard deviations of the differences for GTO minus 
SBUV vary between 0.4 and 0.6%, and the corresponding 
values for ground-based minus SBUV are larger by a factor 
of 1.8 to 2.6. The anomalies inferred from GTO and SBUV 
show very good agreement, with ranges of the differences 
varying between 2.2 and 3.5 %. The corresponding ranges 
for ground-based minus SBUV exhibit larger scattering with 
values ranging between 4.4 and 6.3 %. 

3.4 Seasonal mean zonal mean and anomaly 

It is generally acknowledged that the poor spatial sampling 
and the relatively infrequent measurements could be the ma- 
jor factors causing the differences between ground-based 
and satellite monthly mean zonal mean total ozone data. To 
further explore evidence supporting such an argument, we 
have conducted additional analysis for the comparisons of 
GTO versus SBUV and ground-based versus SBUV using 
the seasonal mean zonal means instead of the monthly mean 
zonal means. For each of the three data records, 61 sea- 
sonal mean zonal means were computed starting from MAM 



1995 2000 2005 2010 

16-year period (MAR 1996 - JUN 201 1| 


GTO minus SBUV Ground-based minus SBUV 

Figure 8. Differences in monthly mean zonal mean anomaly. 

(March- April-May) in 1996 until MAM in 2011. Results of 
the differences in seasonal mean zonal means and in seasonal 
mean zonal mean anomalies are summarized in Tables 6 and 
7, respectively. The entries with brackets in these tables de- 
note the reduction of each parameter compared to the cor- 
responding values in Tables 4 and 5, which are based on 
monthly mean zonal means. The results in Tables 6 and 7 
clearly indicate that the reduction in both the standard devia- 
tions and the ranges of the differences for ground-based mi- 
nus SBUV, when switching from monthly mean zonal means 
to seasonal mean zonal means, are found to significantly ex- 
ceed the corresponding reduction for GTO minus SBUV. 

These results have led to the conclusion that ground-based 
zonal means, while showing no systematic differences with 
respect to satellite-based data records, show larger scatter in 
monthly mean data compared to satellite-based records. The 
differences in the scattering are significantly reduced if sea- 
sonal zonal averages are analyzed. 

3.5 Trend of differences with respect to satellite data 
records 

Another approach to examine the consistency in the variabil- 
ity of total ozone exhibited by the multi-year data records 
is to investigate whether there is any time dependency in 
their differences. In this regard, we have computed the linear 
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Table 5. Summary of the differences in monthly zonal mean anomaly (four broad latitudinal zones). 


Latitudinal 

zones 

GTO minus SBUV (in %) 

Ground-based minus SBUV (in %) 

Standard 

deviation 

Range 

(Max-Min) 

Standard 

deviation 

Range 

(Max-Min) 

(0-30° S) 

0.45 

2.89 

0.79 

4.35 

(0-30° N) 

0.40 

2.19 

1.06 

6.29 

(50-30° S) 

0.59 

3.53 

1.16 

5.74 

(30-60° N) 

0.47 

3.19 

0.87 

5.79 


Table 6. Summary of the differences in seasonal zonal mean total ozone (four broad latitudinal zones). 


Latitudinal GTO minus SBUV (in %) Ground-based minus SBUV (in %) 


zones 

Mean 

Standard 

deviation 

Range 

(Max-Min) 

Mean 

Standard 

deviation 

Range 

(Max-Min) 

(0-30° S) 

0.76 

0.55 

2.72 

0.56 

0.73 

3.82 



(0.08) 

(0.91) 


(0.18) 

(1.32) 

(0-30° N) 

0.43 

0.48 

1.99 

0.10 

0.97 

4.60 



(0.10) 

(0.84) 


(0.15) 

(2.33) 

(50-30° S) 

0.67 

0.49 

2.74 

-0.33 

0.89 

4.51 



(0.17) 

(1.08) 


(0.34) 

(1.25) 

(30-60° N) 

0.20 

0.54 

2.21 

0.10 

0.79 

3.66 



(0.09) 

(1.14) 


(0.19) 

(2.65) 


The entries in brackets represent the reduction compared to analysis based on monthly zonal means. 


trends using the differences of monthly mean zonal means 
for GTO minus SBUV and ground-based minus SBUV rep- 
resented by the red curve and black curve depicted in Fig. 6. 
The trends of differences for the four zonal bands are listed 
in Table 8a in terms of percent per year. The corresponding 
trends in DU per year are listed in Table 8b. These results 
have provided strong evidence that both the differences GTO 
minus SBUV and ground-based minus SBUV show no sig- 
nificant trends for the 16-year period (1996-2011) under in- 
vestigation. 

Quantitatively speaking, the trends are found to vary be- 
tween —0.04 and 0.1 %yr -1 or —0.1 and 0.3DUyr _1 . The 
slightly larger trends for the differences of ground-based mi- 
nus SBUV at 0-30° N are mainly driven by the high bias of 
ground-based data points in early 2010 as indicated by the 
black curve in Fig. 6b. Future updates of the ground-based 
data record might achieve further improvement in this aspect. 

3.6 Deviations from pre-1980 level 

Even though the investigation of ozone trends is beyond the 
scope of this paper, it is still interesting to examine the de- 
viations of total column ozone from pre-1980 levels inferred 
from the three multi-year data records. Since a substantial 
part of ozone variability is related to QBO and the 11 -year 
solar cycle (Bowman, 1989; Hamilton, 1989; Bojkov and 
Fioletov, 1996), removing these natural components from 


the data records will make it easier to examine the long-term 
changes. Our analysis to achieve this purpose was conducted 
as follows: 

i. A regression model fit was performed using the ground- 
based monthly zonal mean total ozone time series cov- 
ering 1964-2011. The model includes the annual cy- 
cle, EESC (effective equivalent stratospheric chlorine)- 
related trend, solar-cycle-related component (using so- 
lar flux at 10.7 cm), QBO-related component (using the 
normalized equatorial wind at 30 and 50 hPa), and vol- 
canic component (using stratospheric optical depth at 
550 nm, only for the year following El Chichon and 
Pinatubo eruptions). 

ii. The QBO-related and solar-cycle-related variations ob- 
tained from step (i) were subtracted from all three data 
records. 

The annual deviations from pre-1980 levels are computed 
and depicted in Fig. 9a-c for 30-60° N, 50-30° S, and 30° S- 
30° N, respectively. Pre-1980 annual zonal means based on 
1964-1980 ground-based data record are used as our base- 
line. 

It is noted that all three data records show consistent pat- 
terns of year-to-year changes for the three zonal bands. 


www.atmos-meas-tech.net/7/1 68 1/2014/ 


Atmos. Meas. Tech., 7, 1681-1692, 2014 



1690 


E. W. Chiou et al.: Comparison of profile total ozone from SBUV (v8.6) 


Table 7. Summary of the differences in seasonal zonal mean anomaly (four broad latitudinal zones). 


Latitudinal 

zones 

GTO minus SBUV (in %) 

Ground-based minus SBUV (in %) 

Standard 

deviation 

Range 

(Max-Min) 

Standard 

deviation 

Range 

(Max-Min) 

(0-30° S) 

0.38 

1.76 

0.62 

3.09 


(0.07) 

(1.13) 

(0.17) 

(1.26) 

(0-30° N) 

0.34 

1.57 

0.92 

4.11 


(0.06) 

(0.62) 

(0.14) 

(2.18) 

(50-30° S) 

0.51 

2.74 

0.89 

4.51 


(0.08) 

(0.79) 

(0.27) 

(1.23) 

(30-60° N) 

0.38 

2.08 

0.71 

3.70 


(0.09) 

(1.11) 

(0.16) 

(2.09) 


The entries in brackets represent the reduction compared to analysis based on monthly zonal mean anomalies. 


Table 8. (a) Trend of the differences (in % per year) for monthly zonal mean total ozone and (b) trend of the differences (in DU per year) for 
monthly zonal mean total ozone. 


(a) 

Latitudinal zone 

GTO minus SBUV 

Ground-based minus SBUV 

(0-30° S) 
(0-30° N) 
(50-30° S) 
(30-60° N) 

0.05 ± 0.0098 % yr _ 1 
0.05 ± 0.0090 % yr _ 1 
0.07 ± 0.0097 %yr _1 
0.04 ± 0.0099 %yr _1 

-0.04 ± 0.0148 %yr _1 
0.12 ± 0.0160 %yr“* 
0.04 ± 0.0204 %yr _l 
-0.02 ± 0.0160 %yr“‘ 

(b) 

Latitudinal zone 

GTO minus SBUV 

Ground-based minus SBUV 

(0-30° S) 
(0-30° N) 
(50-30° S) 
(30-60° N) 

0.14 ±0.0260 DU yr _l 
0.14 ±0.0245 DU yr _ 1 
0.22 ±0.0293 DU yr _ 1 
0.14 ±0.0330 DU yr _l 

-0.11 ±0.0387 DU yr _ 1 
0.31 ±0.0430 DU yr“* 
0.12 ± 0.0617 DU yr -1 
-0.07 ±0.0540 DU yr -1 


For the overlap period of 1996 to 2010, all three data 
records indicate a gradual increase at 30-60° N from —5 % 
in 1996 to —2 % in 2010. In contrast, at 50-30° S and 30° S- 
30° N there has been a levelling off in the 15 years after 
1996. The differences in the deviations inferred from GTO 
and SBUV are depicted in Fig. 10. The results indicate that 
the differences are between 0 % and 1 % with the exceptions 
of the dips in 1996 and 2001. In addition, there is a slight 
increase in the differences during the 15-year period (1996- 
2010 ). 


4 Concluding remarks 

We have presented comparisons of zonal mean total col- 
umn ozone inferred from three recently released indepen- 
dent multi-year data records. The analyses are based on the 
16-year overlap period (March 1996 through June 2011). 
The results of our investigation have led to the conclu- 
sion that, despite the differences in the satellite sensors and 


retrieval methods, the SBUV v8.6 merged profile total ozone 
and GTO merged total ozone data records show very good 
agreement in the monthly mean zonal mean total ozone 
and monthly mean zonal mean anomalies. The ground-based 
zonal means, while showing no systematic differences with 
respect to satellite data records, show larger scatter in the 
monthly mean data compared to satellite-based records. The 
scattering is significantly reduced if seasonal zonal averages 
are analyzed. 

The major findings based on the characteristics in the four 
latitudinal zones 0-30° S, 0-30° N, 50-30° S, and 30-60° N 
can be summarized as follows: 

1. It has been found that, on average, the differences in 
monthly zonal mean total ozone vary between —0.3 and 
0.8 % and are well within 1 %. 

2. The standard deviations of the differences in monthly 
zonal mean total ozone for (GTO minus SBUV) vary 
between 0.6 and 0.7 %. The corresponding differences 
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Deviations from Pre-1980 level 




GTO (GOME Type Total Ozone) * * SBUV(v8.6) 

Ground-based Records 


Figure 9. Deviations in annual mean total ozone from pre-1980 lev- 
els: (a) 30°N-60° N, (b) 50° S-30°S, and (c) 30° S-30° N. 

for ground-based minus SBUV are larger by a factor of 

1.4 to 1.9. 

3. The ranges of the differences GTO minus SBUV vary 
between 2.8 and 3.8 %, while the corresponding ranges 
for ground-based minus SBUV are larger by a factor of 

1.4 to 2.4 with values ranging between 5.1 and 6 . 9 %. 


Differences in the Deviations from Pre-1980 level 



1995 2000 2005 2010 


30N to 60N 50S to 30S 

30S to 30N 

Figure 10. Differences in the deviations from pre-1980 level (GTO- 
SBUV). 

Analyses of the annual deviations from pre-1980 level in- 
dicate that, for the period of 1996 to 2010, all three data 
records show a gradual increase at 30-60° N from —5 % in 
1996 to -2% in 2010. In contrast, at 50-30° S and 30° S- 
30° N there has been a levelling off in the 15 years after 1996. 
The deviations inferred from GTO and SBUV show agree- 
ment within 1 %, but a slight increase has been found in the 
differences during the 15-year period (1996-2010). 
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4. The standard deviations of differences GTO minus 
SBUV in monthly zonal mean anomalies vary between 
0.4 and 0.6 %. The corresponding standard deviations 
for ground-based minus SBUV are larger by a factor of 
1.8 to 2.6. 

5. The ranges of the differences in monthly zonal mean 
anomalies for GTO minus SBUV vary between 2.2 and 

3.5 %, while the corresponding ranges for ground-based 
minus SBUV are larger by a factor of 1.5 to 2.9 with 
values ranging between 4.4 and 6.3 %. 

Both the differences GTO minus SBUV and ground-based 
minus SBUV show no significant trends for the 16-year pe- 
riod indicating the absence of time -dependent differences 
among the three data records. 
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